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1.  Introduction 

OLVER  [12]  has  developed  error  bounds  for  the  so-called 
LIOUVILLE-GREEN1 5 (LG)  approximation  [5,  10]  to  the  solution  of 
the  differential  equation 

- J(t)x.  (1) 

His  results  extended  earlier  work  by  BLUMENTHAL  [1] . The  LG 
approximation  is  of  particular  importance  in  applied  mathematics 
because  it  involves  only  elementary  functions.  As  OLVER  [12] 
has  pointed  out,  however,  the  development  of  strict  upper  bounds 
for  the  errors  in  the  approximate  solutions  has  received  rela- 
tively little  attention.  Moreover,  the  error  bounds  that  have 
been  previously  developed  [12,  13]  are  for  boundary- value  prob- 
lems, since  most  problems  of  interest  in  mathematical  physics 
are  boundary-value  problems.  These  error  bounds  do  not  apply 
to  initial-value  problems. 

Thus,  the  purpose  of  this  paper  is  to  develop  error  bounds 
for  the  LG  approximation  to  initial-value  problems.  Furthermore, 
in  recasting  OLVER' s results  in  a form  suitable  for  initial-value 
problems,  we  have  been  able  to  sharpen  his  bounds  for  boundary- 
value  problems.  Although  we  will  develop  our  results  within 
the  context  of  a specific  problem  in  operations  research,  they 
are  clearly  applicable  to  the  general  second  order  initial-value 
problem,  and  it  is  a straightforward  task  to  so  recast  them. 

^ Also  called  the  WKB  approximation  [13]. 


Probably  the  most  widely  used  (at  least  in  the  United 


States)  deterministic  differential-equation  model  in  operations 

2) 

research  (4,15-17]  are  LANCHESTER-type  equations  of  warfare  , 
which  are  named  for  the  pioneering  1914  work  of  F.  W.  LANCHESTER 
[9].  In  this  paper  we  consider  such  a linear,  variable-coeffi- 
cient differential-equation  model  for  combat  between  two  homogene- 
ous forces.  This  model  yields  an  X force-level  equation  equiva- 
lent to  (1) . Unfortunately,  for  even  the  simplest  time-varying 
attrition-rate  coefficients  of  interest,  the  solutions  cannot 
in  general  be  expressed  in  terms  of  either  elementary  functions 
or  tabulated  higher  transcendental  functions  [17] . It  is  there- 
fore natural  to  seek  an  approximation  in  terms  of  elementary 
functions.  Moreover,  one  is  always  interested  in  a simple 
a priori  estimate  for  the  error  in  the  approximate  solution. 

In  this  paper  we  develop  an  erro^  bound  that  is  both  realistic 
and  easy  to  evaluate.  This  bo;,>-4  Xs  then  computed  for  some 
particular  attrition-rate  coefficients  of  interest. 

2.  Variable-Coefficient  LANCHESTER-Type  Equations  of  Modern 
Warfare 

We  consider  the  following  variable-coefficient  LANCHESTER- 
type  equations  of  modern  warf  for  combat  between  two 

homogeneous  forces 

2) 

We  will  refer  to  any  differential-equation  model  of  combat  as 
being  LANCHESTER-type  equations.  The  state  variables  are 
typically  the  force  levels  of  the  various  weapon  system  types. 

3) 

The  term  "of  modern  warfare"  denotes  that  we  are  considering 
linear  differential  equations.  There  are  other  nonlinear  types 
of  LANCHESTER  equations  [4,16]. 
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dx/dt  * -a(t)  y,  dy/dt  * -b(t)x,  (2) 

with  initial  conditions 

x(t»0)  = Xq  , y(t»0)  = yQ 

where  t = 0 denotes  the  time  at  which  the  battle  begins,  x(t) 
and  y(t)  denote  the  numbers  of  X and  Y at  time  t,  and 
a(t)  and  b(t)  denote  time-dependent  LANCHESTER  attrition- 
rate  coefficients.  These  equations  (2)  have  been  hypothesized 
to  model  combat  in  which  both  sides  use  aimed  fire  and  target 
acquisition  times  are  independent  of  the  numbers  of  firers  and 
targets  [2,191.  The  attrition-rate  coefficients  represent  the 
effectiveness  of  each  side's  fire  (i.e.  its  firepower).  Temporal 
variations  in  a side's  fire  effectiveness  (caused  by,  for  example, 
changes  in  force  separation,  combatant  postures,  target  acquisi- 
tion rates,  etc.)  are  modelled  by  the  time-dependent  attrition- 
rate  coefficients.  Further  discussions  of  the  physical  assump- 
tions hypothesized  to  yield  (2) , estimation  of  the  attrition- 
rate  coefficients,  and  the  importance  of  (2)  to  military  opera- 
tions research  are  found  in  [16,171. 

We  assume  that  a(t)  and  b(t)  are  positive  and  twice 
differentiable  for  tg < t < +»  with  tg  < 0.  We  also  assume 
that  a(t)  , b(t)  eL(tg,T)  for  any  finite  T.  We  further  take 
a(t)  and  b(t)  to  be  given  in  the  form  a(t)  ■ k g(t), 

C I 

b(t)  * k^Mt),  where  kft,  k^  are  positive  constants  chosen  so 
that  a(t)/b(t)  ■ kaA^  when  g(t)  ■ h(t)  for  all  t.  We 
introduce  (see  [171)  the  intensity  of  combat  I(t)  and  the 
relative  fire  effectiveness  R(t)  defined  by 

I(t)  - /a (t) b (t) , and  R(t)  - a(t)/b(t).  (3) 
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xi s and  xr  * W (4) 


A large  class  of  tactical  situations  of  interest  can  be 
modelled  with  the  following  general  power  attrition-rate 
coefficients  [17] 

a(t)  - ka(t+C)M,  and  b(t)  * kb(t+C+A)v'  (5) 

where  A,C  >0.  We  will  call  A the  offset  parameter,  since  it 
allows  us  to  model  (with  y,v>0)  battles  between  weapon  systems 
with  different  maximum  effective  ranges.  We  will  call  C the 
starting  parameter,  since  it  allows  us  to  model  (again,  with 
y,v>0)  battles  that  begin  within  the  minimum  of  the  maximum 
effective  ranges  of  the  two  systems.  The  offset  and  starting 
parameters  are  related  to  various  physical  quantities  in  [17] . 

We  observe  that  tQ  = -C.  Also,  a(t),  b(t)  eL(tQ,T)  implies 
that  y, v > -1. 

Prom  (2)  we  may  obtain  the  X force- level  equation 

lt^’  {iTET  It}  a<t)b(t)x  - 0,  (6) 

with  initial  conditions 

x(t°0)  » xQ , and  { [1/a (t) ]dx/dt }t-0  - -yQ. 

We  may  consider  that  tQ  » max(tQ,tQ),  where  t£  denotes  the 
right-most  finite  singularity  of  the  X force-level  equation. 
Furthermore,  we  set  tg * 0 if  there  are  no  finite  singularities. 
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3.  LIOUVILLE-GREEN  Approximation  to  lANCHESTER-Type  Equations 
of  Modern  Warfare  3 


I (s) ds , 


and  denote  x(t»0)  as  xQ.  Then  tq  - 0 for  to~°*  From 
a(t),  b(t)  e L(tjj#T)  it  follows  that  x ■ x(t)  is  well  defined 
by  the  CAUCHY-SCHWARZ  inequality  for  integrals.  The  transforma- 
tion is  also  easily  seen  to  be  invertable.  We  observe  that 
x - xn  is  related  to  the  average  intensity  of  combat  I (t)  by 


I (s) ds> t * tf (t) . 


We  may  call  x - xQ  "the  elapsed  normalized  battle  time,"  since 
the  trans formation  (7)  reparameterizes  the  battle's  time  scale 
in  terms  of  elapsed  time  of  battle  and  average  combat  intensity 
as  (8)  shows  us.  The  substitution  (7)  transforms  (6)  into 


Hi 


with  initial  conditions 


x(x»x0)  - xQ,  and  { [1/R1/2 (t) J dx/dx}T=T  » -yQ. 


Remark  1:  It  is  easily  shown  (see  [17])  that  (9)  may  be 

transformed  into  a linear  second  order  differential  equation 
with  constant  coefficients  if  and  only  if 

JJtT  * constant. 
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transforms  (9)  into  LIOUVILLE's  normal  form  (see  INCE  [7]  or 
KAMKE  [8]) 


{1+F(t>  }X  = 0f  (11) 

with  initial  conditions 

X(t=t0)  * xQ , and  dX/dx(x*xQ)  - “ xoeo ' 

where 

F (t)  * P"  ( t)  /P  (t)  , P(t)  - [R(t)]“1/4,  (12) 


E(t)  * 4lftT  5T£nR'  (13) 

eQ  denotes  e (t*0) , and  P' (x)  denotes  dP/dx. 

Writing  (11)  as  d2X/dx2  - X = F(x)X,  we  may  use  variation 
of  parameters  to  obtain  the  solution  to  (11)  as 


X(x)  * xQ  cosh  (x-xQ)  - (yQ./RjJ+x0e0)sinh(x-x0) 


+ F (a) sinh  (x-o)X(c)do . (14) 


If  one  drops  the  integral  term  in  (14),  one  obtains  the 


LIOUVILLE-GREEN  approximation  X(x) 

4) 

Heuristically , if  the  appropriate  fractional  power  of  the  relative 
effectiveness  R(t)  is  "slowly  varying,"  then  by  (12)  we  would 
expect  that  | F ( x) | <<1  so  that  the  integral  term  in  (14)  is 
"negligible."  Theorem  1 gives  us  bounds  on  how  "negligible" 
this  term  is. 
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t 


(15) 


&(t)  = xQ  cosh  (t-t0)  - (y0/^+x0e0)sinh(T-T0) 


which  in  terms  of  the  original  independent  variable  x reads 
x(t)  - {xQcosh  (t-t0)  - (y0/R^+x0e0)sinh(T-T0)}.  (15) 


We  observe  from  (14)  that  F(t)  >0  for  all  T - tq  implies  that 
as  long  as  x(t)  >0  we  have  x(t)  >x(t).  A similar  statement 
holds  for  F (t)  <0.  As  we  shall  see  below,  such  cases  in  which 
F(t)  is  always  >0  or  £ 0 are  readily  encountered  in  applica- 
tions . 


4 . Error  Bounds  for  the  LI OUVILLE- GREEN  Approximation 
The  main  result  of  this  paper  is  Theorem  1. 


Theorem  1:  Error  bounds  for  the  LI OUVILLE -GREEN  approximation 
are  given  by 

|x(t)-x(t)  | < xQ  Kj  e(t)  < Xg  Ky  e (t)  , (17) 

where 

y 

Ky  - 2 { (It  | Eg  |)  ^ , (18) 

y y 

J»I  for  l-^^SEg  and  Kj  - 1 + Eg  + ^ ^ , (19) 

y y 

J » II  for  "I  - ~ ^ ^ < 1-^  ^ and  KJZ  » 2,  (20) 

y y 

J - HI  for  e0<-1-^^  and  “m  * 1 “ e0  “ 3^  *^0  * 0 ' (21) 

and 
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e(t)  = [^^-]1/4{exp(j|X  |F(a)  |do)-l}sinh(x-T0)  . (22) 

0 ^ 

The  sign  of  the  error  is  determined  by  the  sign  of  F(t) . As 
long  as  x(t)  >0,  it  follows  that 

F(t)  £ 0 for  all  t£0  implies  that  x(t)>x(t), 

with  the  last  inequality  being  reversed  when  F(t)  *0. 

Proof:  Theorem  1 readily  follows  from  Lemma  2,  which  is  proven  below. 

5 . Development  of  Error  Bounds 

Consider  the  following  fundamental  system  of  solutions 
{X1,X2>  to  (11) 

Xk(x)  * {l+hk(T)}exp[(-l)(k"1)  (t-t0)]  for  k - 1,2,  (23> 

where  h^fx)  for  k = 1,2  is  to  be  chosen  so  that 

hk(T=*T0)  = 0#  and  dh^/dxfx^Xg)  =»  0.  (24) 

It  follows  that  the  solution  to  (11)  may  be  expressed  as 
1 <t-t0) 

X(x)  = j{(x0(l-e0)-y0/R^)e  Cl+hj^ (t)  3 

- (X-X  ) 

+ (x0(,l+e0)+y0/R^)e  0 (l+h2(r)]},  (25) 

so  that 

* i (t-tJ 

X(t)-X(t)  = j{  [xQ  (l-e0)-y0/R^Je  uh1(x) 

- (x-x0) 

+ Ix0(l+e0)+y0/R^]e  u h2(x)}.  (26) 

# 
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Substituting  (23)  into  (11),  we  find  that  h^ft)  for  k = 1,2 
roust  satisfy 

d2h  k able 

-jY#--  <-l>*  2-af--F<T)h|£  - F<t>,  (27) 

with  initial  conditions  (24) . 

We  may  consider  hk<T)  for  k - 1,2  to  be  an  error  term 
for  the  LIOUVILLE-GREEN  approximation.  We  next  develop  a bound 
on  its  magnitude,  which  sharpens  earlier  results  by  OLVER  [12]. 


Lemma  Is  A bound  on  the  magnitude  of  hk(T)  for  k = 1,2  is 

given  .Jay 

|hk(T)  |exp  [(-l)(k“1)(T-T0)]  < 2{exp(i|X  |P (a)  |da) -l}sinh (t-tq)  . (28) 

T0 


Proof:  For  notational  convenience,  we  develop  the  bounds  for  (t)  and 
h2 (t)  separately.  Transposing  the  right-most  term  on  the  left-hand  side  of 
(27)  for  k ■ 1,  treating  (l+h^ (t) )f(t)  as  a "forcing  term,"  and  inte- 
grating twice;  we  obtain  the  following  VOLTERKA  integral  equation  after  a 
further  integration  by  parts 


'l(T)  "if  t1“«2(0"T)>f<a>^+h1(a)}do. 

J j 


(29) 


Solving  (29)  in  the  usual  manner  by  successive  approximations,  we  obtain 


h (,)  - I T (,), 
n«l 


(30) 


where  Tq(t) 


1 and  for  nil 


P(t)  - U-e2(°"T)}F(o)T  . (o ) do . 

n 2Jx  n-1 


(31) 
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J 


2<  j -2(t-t  ) 

Observing  that  { 1-e  ' } < l - e 


for  0 < Tq  — o — x , we  find  that 


1ti(t)  I < if  |p(o)|do, 

T0 


with  equality  holding  for  t ■ A straightforward  inductive  argument 

along  the  usual  lines  now  yields 


“2<T-V}  ire  ) n 

|TjT)| “jj^  |F(a)|do|  . 


_n 

2 nl 


(32) 


One  step  in  the  inductive  proof  of  (32)  is  deserving  of  further  elaboration, 
however.  From  (31)  and  (32)  we  obtain 

lTntl<T)l5~ntl~  ~ } U-e2(o‘T>}{l-e  <0  0>}  |f(o)  |{£  |f(u)  |du}  do.  (33) 

The  inductive  proof  of  (32)  is  coopleted  by  combining  (33)  with  the  observa- 

/ -2(o-t  » -2(t-t  ) 

tion  that  { 1-e  ' } ^ 1-e  / 5 1 - e for  0 < tq  5 a 5 t . Our 

sharpening  of  OLVER's  results  is  due  to  this  observation.  The  remaining 
steps  in  the  proof  of  (28)  for  k - 1 follow  along  well-known  lines  [12,13] 
and  will  be  omitted  here.  Similar  arguments  are  used  to  prove  (28)  for 
k - 2. 


5) 


Remark  2:  In  our  notation,  OLVER's  [12]  corresponding 
error  bound  for  k - 1 would  read 

, t 

|h  (t)  | < exp (7/  |f(o)  |do)-l  . 

1 “ * T0 

His  corresponding  error  bound  for  k>2  is  not  directly 
comparable  to  our  result  here,  since  he  does  not  take 
both  errors  zero  at  the  same  point. ^ 


(34) 


In  our  notation  for  the  finite  or  infinite  interval  ^To,Ti*' 
OLVER  [12,13]  takes  h^  (t«t Q)  - d^/difT-t^  - 0 and 
^2^t=t1^  * dhj/dt (t=t^)  » 0 in  contrast  to  (24).  This  is 
what  makes  his  results  unsuitable  for  initial-value  problems. 
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Lemma  2:  In  terms  of  the  transformed  dependent  variable  X(t) , 


error  bounds  for  the  LIOUVILLE-GREEN  approximation  are  given  bj 


|X(t)-X(t)  | < xQ  KjE(t)  < Xq  Ky  E (t)  , 


where  Ky  and  Kj  are  given  by  (18)  through  (21) , and 


E(x)  = (exp( 


\\\  |F<°>I 

l A 


do)-l>sinh (t-Tq) . 


Proof:  From  (26)  and  Lemma  1 we  obtain 


|x(t)-x(t)  | < {|x0(l-eQ)-y0^|  + lx0(l+e0)+y0^|  }e(t). 


It  follows  that  a rather  loose  error  bound  is  given  by 


|x(t)-X(t) | S xQ  Ky  E (t) . (38) 

y 0/—  y0  /— 

We  observe  that  for  -1 /R_  < e_  < 1- — /R_,  we  have 

x0  “ 0 xo  0 

lxo(1*eo)"yo^  + lxo(1+Eo)+yo^  " 2x0'  (39) 

so  that  (37)  becomes  |x(t)-X(t) | 5 2xqE(t) . Thus,  (35)  is  proven  for  J ■ II. 

y0  r- 

For  1-  — /R  <e_,  the  error  bound  (37)  becomes 
x0  o o 

a y 

|x(t)-X(t)|  <2xQ(e0+j2.  /r^)e(t),  (40) 

since 


IVW-yo^l  + lxo(1+co)+yo^  * 2xo(co  + ^ ^>* 
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The  error  bound  (40)  may  be  sharpened,  however,  as  follows.  If  F(o)i0  for 
0<  Tq  — O — t,  then  as  long  as  X(x)>0  we  have  X(x)>X(x)  and  gk(x)>0, 
where  for  k * 1,2 


gk<T)  - jhfc(x)  expI(-l)k(T-T0)l. 


It  follows  from  (26)  that  for  1 — - /r~ < t 

xQ  0 0 


(42) 


where 


0<X(T)  -x(t)  <{x0(l+cQ)  +y0»^>E+(T), 


e+(t)  - {exp(^  | F(a)do)  - l}sinh(T-r0) 


(43) 


(44) 


Similarly,  if  F(o)2  0 for  Tflia<T,  then  as  long  as  X(x)>0  we  have 
X(x)  £X(x)  and  g^(x)  SO,  so  that 


o>x(x)  -x(x)  > (x0n-+€0)+y0’^>E_(x) , 


where 


E_ (x)  » {1-exp ( 


F(o)d a) }sinh(x-xQ) 


(45) 


(46) 


It  follows  from  (43)  and  (4S)  that 


|x(x)  - X(x)  1 2 xn{l  + cn  + ~ i^)E(x). 

0 


(47) 


Jr  y y 

Furthermore,  since  e > 1 - -2-  *''5“  implies  that  2(en  + — > 1 + e_+  — /r~. 


'0--  xQ  "0 


■0  xQ  "O' 


-0  xQ  "O' 


the  bound  given  by  (47)  is  sharper  than  that  given  by  (40).  Thus,  (35)  is 

proven  for  J - I.  The  proof  of  (35)  for  J - III  is  similar  to  that  for 

J ■ II.  It  is  easily  seen  that  K, < K for  x.,y  > 0. 

J u 0 0 


1 


1 

I 


1 

i 


12 


6 . Examples 


We  now  compute  theoretical  error  bounds  for  two  special 
cases  of  the  LIOUVILLE-GREEN  approximation  to  the  solution  of 
(6)  with  the  general  power  attrition-rate  coefficients  (5) : 

(I)  power  attrition-rate  coefficients  with  no  offset  (i.e.  A»0), 

and  (II)  linear  attrition-rate  coefficients  with  positive  offset 
(i.e.  A > 0)  . 

6.1.  Power  Attrition-Rate  Coefficients  with  No  Offset 
In  this  case  we  have 

a(t)  * ka(t+C)y,  and  b(t)  * kb(t+C)v,  (48) 

with  C>0  and  y,v > -1.  The  LIOUVILLE-GREEN  approximation  to 
the  X force  level  is 

x ( t)  » (1+t/C)  *y_v^4{x0  cosh  (t-Tq) 

c(M-v)/2  + _0_^ c“5]sinh(T-T0)  } , (49) 

where 

x(t)  - (l/fiJXjtt+C)6,  (50) 

and 

6 - (m+v+2)/2.  (51) 

In  preparation  for  estimating  the  error  in  the  LIOUVILLE- 
GREEN  approximation  (49)  by  Theorem  1,  we  compute 

P(t)  - iHiy)  < . (32' 

166*t* 
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Thus,  F(t)  >0  for  all  t > > 0 if  and  only  if  \i  2 v.  Hence, 

as  noted  at  the  end  of  Section  3,  one  frequently  encounters  in 
applications  cases  in  which  F(t)  always  has  the  same  sign. 

For  the  error  estimate  (17) , we  then  have 

| |T  |F(o)  |do  - ll!^^±i±iL{c"6-(t+C)'6}.  (53) 

Remark  4:  The  exact  solution  x(t)  to  (6)  with  attrition- 

rate  coefficients  (48)  is  given  in  TAYLOR  and  BROWN  [17]  (see 

also  [16]).  It  may  be  expressed  in  terms  of  modified  Bessel 

functions  of  the  first  kind  of  (for  p ,v  > -1)  fractional  order, 

i.e.  I for  0 < a < 1.  Since  few  of  the  latter  are  tabulated 
Cl 

(i.e.  tabulations  only  exist  for  a ■ ±1/4,  ±1/3,  ±1/2,  ±2/3, 

± 3/4,  and  these  do  correspond  to  cases  of  interest) , TAYLOR 

and  BROWN  [17]  suggested  the  use  of  new  transcendents  which 

they  called  LANCHESTER-CLIFFORD-SCHLAFLI  functions6* . 

?• 

J 

6.2.  Linear  Attrition-Rate  Coefficients  with  Positive 
Offset 

a (t)  * ka(t+C),  and  b(t)  * kfa(t+C+A),  (54) 

with  A,C>0.  The  LIOUVILLE-GREEN  approximation  to  the  X force 
level  is 


6) 


After  earlier  work  by  W.  K.  CLIFFORD  [3]  and  L.  SCHLAFLI  [14) 
(see  also  [6,18] ) . 


x(t)  ■ [(llv^c))]1/<txoco8h  (T-To) 


. ~yo/T»i  xoA/g 

L/l+A/C  4XjC2 (l+A/C) 3/2 


™sinh(x-x0)}. 


(55) 


where 


x (t)  - Jln(<|>+/iF^I)}, 


(56) 


and 


<Mt)  - 1 + 2(t+C)/A. 


(57) 


For  estimating  the  error  in  the  LIOUVILLE-GREEN  approxi- 
mation (55),  it  is  more  convenient  to  express  F(t)  as  defined 


by  (12)  in  terms  of  the  original  independent  variable  t.  We 

.7) 


compute  that 


_ A(l2(t+C)  + 7A> 

F(x)  - n * ? . 

161^ ( t+C)  (t+C+A) J 


(58) 


It  follows  that  F(t)  > 0 for  all  x » Tq  > 0.  Thus,  as  in  the 


previous  example,  F(x)  always  has  the  same  sign.  For  the 
error  estimate  (17) , we  then  have 


F (o)  do  < min  (n^  (t)  ,m2  (t) ) , 


(59) 


where 


m,  (t) 


3A 


16C  (1+A/C)  ' A 


T7T. 


(4  (l-q1/2)  + ^(l-q3/2)  ) , (60) 


m2  (t) 


32CTA 


3—  { 4 (1-q3)  +2^(l-q4)}. 


(61) 


7) 


In  general 
F (x ) 


^7  A {- 
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and 


q(t)  * l/(l+t/C). 


It  may  be  shown  that  m^(t^)  * m2(t^)  implies  that  m1(t)  >m2(t) 
for  all  t>t..  The  error  control  term  estimate  m^(t)  was 
developed  for  "small"  t8^  , while  m2  (t)  for  "lar9e"  . 

Remark  5 : This  case  is  of  particular  interest  in  military 
operations  research,  since  it  may  be  used  to  study  combat 
between  two  weapon  systems  with  different  maximum  effective 
ranges  [16).  The  exact  solution  x(t)  to  (6)  with  attrition- 
rate  coefficients  (54)  is  given  in  [17]  (see  also  [16]).  It 
apparently  cannot  be  expressed  in  terms  of  previously  "known" 
transcendents,  since  the  X force- level  equation  in  this  case 
could  not  be  found  to  correspond  to  any  second  order  linear 
equation  considered  in  [8]  or  [11]. 


7 . Final  Remarks 

Although  given  within  the  context  of  a specific  problem 
in  operations  research,  the  reader  will  have  no  trouble  trans- 
lating the  above  results  into  those  for  the  general  second  order 


Letting  J 


r {12(s+C)+7A}ds 
0 (s+C)5/2(s+C+A) 5/ 


we  find  that 


■ 71  ■ . Here  we  have  used  the  bound  J s 

16AI  (C+A) 


bn  L 


F (o ) do 


{12  (s+C)+7A}d s 
■ 

(s+C)-' 


In  this  case  (to  be  contrasted  with  the  previous  one),  we  have 


used  the  bound  J 


< {12(s+C)+7A}ds 

~ Jo  (s+C)5 
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initial-value  problem.  In  the  examples  of  Section  6 we  saw 
that  for  two  models  of  considerable  importance  in  military  opera- 
tions research  both  the  LIOUVILLE-GREEN  approximation  and  bounds 
on  its  error  were  simply  expressed  in  terms  of  elementary  func- 
tions. No  previous  application  of  the  LIOUVILLE-GREEN  approxi- 
mation has  appeared  in  the  operations  research  literature.  In 
a subsequent  paper  we  plan  to  present  a numerical  investigation 
of  the  accuracy  (both  numerical  eveluation  of  the  theoretical 
error  bounds  and  a comparison  with  the  exact  solution  when 
available)  of  the  LIOUVILLE-GREEN  approximation  to  the  solution 
of  LANCHESTER-type  equations  of  modern  warfare  for  combat  between 
two  homogeneous  forces. 
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